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ABSTRACT 

A  Z -transform  algorithm,  developed  for  the  spectral  analysis  of  signals,  allows 
one  to  get  closer  to  the  poles  of  a  signal  and  effectively  reduces  the  signal’s  band¬ 
width  and  sharpens  its  peak  point.  It  can  give  a  high  resolution,  narrow-band  fre¬ 
quency  analysis  with  frequency  spacing  Af<.l/T,  where  T  =  total  length  of  the 
analysis  interval.  This  algorithm  also  enhances  (1)  the  signal  poles  that  lie  on  cir¬ 
cular  or  spiral  contours  that  begin  at  almost  any  point  in  the  Z-plane  and  (2)  the 
angular  spacing  of  points  in  an  arbitrary  constant.  Since  this  algorithm  takes  advan¬ 
tage  of  high-speed  convolution,  it  is  almost  as  fast  as  and  more  flexible  than  the 
Fast  Fourier  Transform  (FFT). 
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SPECTRAL  ANALYSIS  OF  SIGNALS  BY  USING 
THE  Z-TRANSFORM  ALGORITHM 


INTRODUCTION 

A  Z-transform  algorithm  is  applied  to  the  spectral  analysis  of 
signals  whose  poles  and  zeros  are  either  on  or  off  the  imaginary  axis 
of  a  complex  S-plane.  This  algorithm  is  flexible  and  can  be  applied 
either  on  or  off  the  imaginary  axis  and  can  enhance  the  poles  in  the 
spectral  analysis  of  signals.  It  is  also  useful  in  evaluating  high  reso¬ 
lution,  narrow -frequency-band  spectra. 

The  signals  that  are  to  be  analyzed  by  using  the  digital  technique  are 
considered  at  discrete  values  of  t,  usually  nAt,  where  n  =  0,  1,  2,  .  .  ., 
and  At  is  a  fixed  positive  number  usually  referred  to  as  a  sampling 
period.  In  Fig.  1,  a  continuous  function  of  time,  f(t),  is  shown,  and 
its  values  at  t  =  nAt  are  indicated. 

Let  At  be  the  fixed  positive  number  and  let  f(t)  be  defined  for 
t  >  0.  The  Z-transform  of  f(t)  is  the  function 

oo 

F (z)  =  2  f(nAt)  Z"n  ,  (1) 

n  =0 

where 

Z  =  complex  variable. 

For  finite  sum,  we  can  write 

N-  1 

F (z)  =  2  f (nAt)  Z~n  .  (2) 

n  =0 

If  Z  =  e  ,  where  S  =  a  +  jto,  then  Eq.  (2)  stands  for  the  Laplace 
transform  of  a  train  of  equally  spaced  impulses  of  magnitude  f(nAt). 
The  Laplace  transform  of  a  train  of  impulses  repeats  its  value  in  a 
horizontal  strip  of  the  S-plane  of  width  Zvc  /  At  in  every  strip  parallel 
to  it.  The  Z-transform  maps  each  such  strip  into  an  entire  Z-plane  or, 
conversely,  the  entire  Z-plane  corresponds  to  any  horizontal  strip 
of  the  S-plane,  e.  g.  ,  the  regions  -  ~  <  a  <  <*>  and  -tt /At  _<  to  <,77 /At, 
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where  S  =  cr  +  jco.  The  jco  axis  of  the  S-plane  along  which  we  generally 
equate  the  Laplace  transform  with  the  Fourier  transform  is  the  unit 
circle  in  the  Z-plane,  and  the  origin  of  the  S-plane  corresponds  to 
Z  =  1.  The  interior  of  the  unit  circle  of  a  Z-plane  corresponds  to  the 
left  half  of  the  S-plane,  and  the  external  corresponds  to  the  right-half 
plane.  Straight  lines  in  the  S-plane  correspond  to  a  circle  or  a  spiral 
in  the  Z-plane,  as  shown  in  Fig.  2. 

Usually  the  Z-transforms  are  computed  along  the  path  corresponding 
to  the  jco  axis;  i.  e.  ,  the  unit  circle  on  the  Z-plane  gives  a  discrete 
equivalent  of  the  Fourier  transform. 

There  are  many  applications  of  the  Z -transform  on  and  off  the 
unit  circle,1’2  e.  g.  ,  spectra,  filtering,  interpolation,  etc.  This  report 
presents  the  spectral  estimation  of  signals  for  which  one  has  some 
a  priori  knowledge  of  the  locations  of  poles  by  using  the  Z-transform 
algorithm.  Moreover,  this  report  investigates  the  effect  of  additive 
white  noise  on  the  Z-transform  algorithm  when  signal  spectrums  are 
evaluated. 


THEORETICAL  APPROACH 

For  convenience,  the  equation  for  the  finite  sum(Eq.  (2))  is  written 
again1,2  : 


n  - 1 

F(zk)  =  2 

n  =  0 


f(nAt)  Zk-n 


(3) 


Let 


Z.  =  AW"k,  k  =0, 1,2,.  •  M-l  , 


(4) 


where  M  is  an  arbitrary  integer  and  both  A  and  W  are  arbitrary 
complex  numbers: 

A  =  AQ  e,2ndQ  (5) 

and 

W=Woei2?70o  .  (6) 

The  starting  point  of  the  contour  is  determined  by  A,  and  the  reso¬ 
lution  is  determined  by  W  (Fig.  3). 


If  A  =  1,  W  =e-i  277/N  ,  and  N  =  M,  then  Eq.  (3)  corresponds  to 
the  discrete  Fourier  transform.  If  N  is  a  power  of  two,  the  algorithm 
can  be  implemented  as  a  Fast  Fourier  Transform  (FFT). 
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The  general  Z  -plane  contour  begins  with  point  Z  =  A  and,  depending 
on  the  value  of  W,  spirals  in  or  out  withrespect  to  its  origin.  If  W0  =  1, 
the  contour  is  an  arc  of  a  circle.  The  angular  spacing  of  the  samples 
in  the  Z -plane  is  2tt4>0. 


Since  A  and  W  are  arbitrary  complex  numbers,  we  see  that  the 
points 

sk  =s0  +  k  (Aa  +  jAco)  =^-(ia  A -kin  W)  ,  (7) 

where 

k  =  0,  1,  2, .  .  .  ,  M  —  1  and 
so  =ao  +i<uo  =  \7ln  A  ’ 

lie  on  an  arbitrary  straight-line  segment  of  arbitrary  length  and  arbitrary 
sampling  density. 

The  contour  =  AW_k,  where  k  =  0,  1,  3,  .  .  .  ,  M-l,  is  con¬ 
siderably  more  general  than  that  to  which  the  FFT  applies  (Fig.  3).  It 
is  seen  from  Eq.  (3)  that  the  computation  of  a  Z-transform  along  the 
general  contour  requires  NM  multiplication  and  additions,  since  the 
special  symmetry  of  e)277k/N  ,  which  is  exploited  in  deriving  the  FFT,  is 
absent  in  more  general  cases.  However,  by  manipulating  the  role  of 
Wn  /2  in  calculation,  it  can  be  shown  that  the  computations  can  be  reduced. 


Substituting  the  value  of  Zfc  =  AW'k  ,  we  see  that  Eq. 

F(Z.)=  £  f(nAt)  (AW-k)“n 

n  =0 
N-  1 

=  2  fn  A~n  Wnk  , 

i  n  =0 

where 

f  n  =  f  (n  A  t)  , 


(3)  yields 

(8) 


n  =  0,  1,  2, .  .  . ,  N  — 1,  and 
k  =0,  1,2,.  .  . ,  M-l  • 


Substituting  the  value  of  nk  =  n2  +  k2  -  (k  -  n)2  /Z  for  the  exponent  of 
W  in  Eq.  (8)  we  get 


F(Zk) 


N-l 

2 

n  =0 


fn  A“n  W 


n2  +  k2 


—  (k  -n)  2 
2 
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=  2  fn  A”"  W”^2  Wk  12  W“(k  n 
n=0 

=  2  fn  A'n  W1^2  (w"(k-n)2//2)  (wk2^2)  . 

n  =  0 

Let 

Yn  =  fn  A~n  W”2^2  ,  n  =  0,  1,  2, .  .  . ,  N  -1  ,  (10) 

and 

Vn  =W'n2//2  ,n  =0,  1,  2,  .  .  .,N-1  .  (11) 

Now  convolving  Yn  with  Vn  ,  we  can  write 

N  -  1 

gk=  S  YnV  k=0,l,2,...,M-l  .  (12) 

n  =  0 

k2AI 

Multiplying  Eq.  (12)  by  W  ,  we  can  write 

F(Zk)  =gk  Wk2/2  ,k  =  0,  1,2,.  .  .,M-1  ,  (13) 

which  is  evaluated  as  follows: 

12/1 

a.  The  input  signals  in  sequence  are  multiplied  by  A~nWn  2  to 
get  Yn  ,  which  requires  N  multiplication. 

b.  The  symbol  gk  is  the  convolution  of  Yn  and  Vn  .  This  con¬ 
volution  is  equivalent  to  multiplying  the  Fourier  transform  of  Yn  and 
Vn  .  The  gk  is  evaluated  by  taking  the  inverse  FFT  of  the  product  of 
the  individual  FFT  of  Yn  and  Vn  .  These  operations  take  approximately 
a  time  proportional  to  (N  +  M)  log  (N  +  M).2 

k  2/2 

c.  The  next  step  is  to  multiply  gk  by  W  to  get  F(Zk),  which 
requires  only  M  multiplication.  (The  whole  operation  is  summarized  in 

Fig.  4.) 


SPECTRAL  ANALYSIS  OF  SIGNALS 

The  advantage  of  the  Z -transform  algorithm  over  the  FFT  is  the 
algorithm's  ability  to  evaluate  the  Z-transform  at  points  inside  and 
outside  the  unit  circle.  By  evaluating  the  transform  that  is  outside  the 
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unit  circle,  one  can  get  closer  to  the  poles  and  the  zeros  of  the  signal, 
thus  effectively  reducing  the  signal's  bandwidth  and  sharpening  its  peak 
point.  This  will  be  demonstrated  by  the  simulated  examples  1  and  2. 

Example  1: 

4 

f(t)=  2  A.  e  a‘  cos  Inl.t,  0  <  t  <  T  =  78  msec, 
i=  1  1 

=  0  T  <t  .<  128  msec, 


where 

A.  =  1,  i  =  1, 

2,  3,  and  4 

cq  =-0.0038 

ft  =  50  Hz 

a2  =  -0.0164 

f2  =  200  Hz 

a3  =  -0.023 

f?  =  300  Hz 

a4  =  -0.029 

f4  =  370  Hz 

and  l/At  =  sampling  frequency  =  1  kHz. 

The  signal  f(t)  is  shown  in  Fig.  5.  For  this  example,  only  the  pole 
positions  of  the  upper  half  of  the  S-plane  are  shown  in  Fig.  6. 

In  Example  1,  the  Z-transform  of  128  sample  datapoints  is  evaluated 
on  two  spirals  outside  the  unit  circle;  one  spiral  is  on  the  unit  circle 
and  the  other  is  inside  the  unit  circle.  Figure  7  shows  the  four  contours 
as  they  would  appear  in  the  S-plane.  The  Z-transform  is  evaluated  at 
51  equally  spaced  points,  from  0  through  500  Hz  (frequency  spacing 
Af  =  10  Hz),  that  have  a  corresponding  value  of  <J>0  =  -l/l00.  The 
evaluation  of  the  magnitude  response  of  simulated  signals  without  noise 
on  contours  1,  2,  3,  and  4  of  Fig.  7  is  shown  in  Figs.  8A  through  8D, 
respectively.  Figures  8A  and8B  show  that  peak  values  at  higher  domi¬ 
nant  frequencies  are  smaller  than  the  peak  values  at  lower  dominant 
frequencies  because  the  signal  poles  are  farther  away  at  higher  fre¬ 
quencies  than  at  lower  frequencies.  Figure  8C  shows  that  the  peak 
values  at  all  dominant  frequencies  are  almost  equal.  This  is  due  to 
the  fact  that  the  gradient  of  contour  4  of  Fig.  7  relative  to  contour  3  is 
constant.  As  expected,  Fig.  8D  shows  that  the  sharpening  of  magnitude 
response  in  the  region  of  the  poles  is  quite  pronounced  and  higher  com¬ 
pared  with  Figs.  8A ,  8B,  and  8C  since  the  poles  are  very  close  to 

contour  4.  For  comparison,  a  direct  FFT  of  simulated  signals  is 
evaluated  and  shown  in  Fig.  9.  The  resolution  or  frequency  spacing  is 
Af  =  l/NAt  =  7.  8  Hz,  whereas  the  frequency  spacing  in  the  case  of  a 
Z-transform  is  10  Hz. 
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To  see  the  effect  of  noise  on  the  Z -transform  algorithm,  white  noise 
is  added  to  the  signal.  Now,  the  Z-transform  of  the  composite  signal 
(signal  plus  noise)  is  evaluated,  as  mentioned  above,  on  the  same  con¬ 
tours  at  signal-to-noise  ratio  (S/N)  =  +8.87  dB  and  at  S/N  =  -0.6dB. 
The  S/N  is  defined  as 

S  j  average  signal  power 

N  average  noise  power 

Figures  10A  through  10D  and  1 1A  through  1  ID,  which  show  the  magni¬ 
tude  response  of  composite  signals  at  S/N  =  8.  87  dB  and  -0.  6  dB, 
respectively,  indicate  that  the  magnitude  response  in  the  pole  regions 
is  more  pronounced.  Since  the  noise  is  white,  its  effect  is  noted  by  the 
bias  increase  of  the  whole  spectrum.  For  comparison,  direct  FFT 
evaluation  of  composite  signals  is  shown  in  Figs.  12A  and  12B. 

Example  2:  One  very  useful  advantage  of  the  Z -transform  algo¬ 
rithm  is  that  it  enables  evaluating  high  resolution,  narrow  frequency 
band  spectra  efficiently,  whereas  the  frequency  resolution  in  the  stand¬ 
ard  FFT  technique  is  limited  by  the  total  length  of  the  signal.  However, 
if  one  requires  resolution,  Af,  where  Af  l/NAt  and  where  At  = 
sampling  period  and  N  =  number  of  sample  points ,  then  the  number  of 
sample  points,  N,  must  be  increased.  In  other  words,  the  small  Af 
implies  a  large  value  of  N. 

Often,  however,  a  high  resolution  for  a  limited  range  of  frequencies 
and  low  resolution  for  the  rest  of  the  spectrum  is  desired,  i.  e.  ,  a 
detailed  view  of  the  spectrum  in  the  passband  and  only  a  general  view 
outside  the  passband. 

The  following  simulated  results  will  show  that  the  Z-transform  algo¬ 
rithm  allows  (1)  independent  selection  of  initial  frequency  or  starting 
point  and  (2)  frequency  spacing  that  is  independent  of  the  number  of 
time  samples.  (For  the  standard  FFT,  the  frequency  spacing  is  not 
independent  of  the  number  of  time  samples.)  Hence,  high  resolution 
can  be  obtained  over  a  narrow  frequency  range. 

To  illustrate  these  points,  a  signal  in  the  form  of 
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f(t)  =  2  cos  277" f  .t,  0  <t  <T  =  78msec, 
i  =  1  1 

=  0  T  <  t  <  128  msec  , 
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where 

fx  =  230  Hz 

f2  =  240  Hz 

f3  =  250  Hz 

f„  =  260  Hz 

f5  =  270  Hz  , 

is  generated  as  shown  in  Fig.  13. 

The  FFT  for  128  points  of  this  signal  is  calculated  at  various  S/N1  s  , 
as  shown  in  Figs.  14A,  14B,  and  14C.  In  all  cases,  most  of  the  signal 
energy  is  located  between  200  through  300  Hz. 

To  investigate  the  passband  and  the  transition  regions  more  care¬ 
fully,  the  Z-transform  algorithm  was  used  to  give  a  2-Hz  resolution 
from  200  through  300  Hz,  as  shown  in  Figs.  15A,  15B,  and  15C.  The 
contour  is  the  same  as  in  the  case  of  the  FFT.  To  achieve  this  reso¬ 
lution,  one  would  require  500  time-sample  values  instead  of  only  the 
128  points  that  have  been  used  in  the  Z-transform  algorithm.  This 
shows  that  this  algorithm  is  a  powerful  tool  for  the  close  examination 
of  small  frequency  bands. 


CONCLUSION 

Signal  simulation  results  show  that  the  Z-transform  algorithm  can 
enhance  the  signal  poles,  which,  in  turn,  sharpen  the  peak  of  the  domi¬ 
nant  frequency  and  reduce  the  bandwidth.  The  algorithm  provides  a 
number  of  output  points  that  are  independent  of  the  number  of  input 
points;  i.  e.  ,  the  output  points  may  be  larger  or  smaller  than  the  input 
points.  It  also  has  the  added  advantage  of  enabling  one  to  start  at  an 
arbitrary  point  in  the  Z-plane.  This  algorithm  will  be  very  useful  for 
a  close  look  at  the  filter  passbands. 
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Fig.  2.  Correspondence  of  a  Z-Plane  Contour  with  an  S-Plane  Contour 


8 


A'nWn2/*  ,n  =  0,l,2,...,N-l 


©-^ 


f  (At)  ,n  =  0,1 1 


HIGH-SPEED 
CONVOLUTION 
Yn  *  w“n2/2 


Ok 


W 


k2/2 


■0- 


F(Zk)  ,k  =  0,l,2,...,M-l 


Fig.  4.  Steps  Involved  in  Evaluating  a  Z-Transform  of  Time  Function 
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f  (t)  =  z  A.eai  cos  2 7T f  j t  ,  0  <t<T  =  78  msec 
i=l 

=  0  T  <t<  128msec 


where  a,=  -0.0038  f,=  50Hz 
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Fig.  5.  Simulated  Signal 
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Fig.  6. 


Pole  Positions  of  a  Simulated  Signal  in  the  Left  Upper-Half  of  an  S-Plane 


11 


12 


wnyiD3dS  do  amvA  aimosav 


00* 

oi) 

E 


i 


13 


Z-Transform  Absolute  Magnitude  Response  of  Simulated  Signals  Evaluated  on  the 
Four  Contours  of  Fig.  7  at  S/N  =oo  dB 
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Fig.  10.  Z-Transform  Absolute  Magnitude  Response  of  Simulated  Signals  Evaluated  on  the 
Four  Contours  of  Fig.  7  at  S/N  =8.87  dB 
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Fig.  11.  Z-Transform  Absolute  Magnitude  Response  of  Simulated  Signals  Evaluated  on  the 
Four  Contours  of  Fig.  7  at  S/N  =—0.6  dB 
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Fig.  12.  FFT  Absolute  Magnitude  of  Simulated  Signals  at  S/N  =8.87  dB  and  S/N  =  —0.6  dB, 
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f  (t)  =  £  cos  2 *7T f j f  #  0  <  t <T  =  78  msec  1 

i  =  t 

=  0  T  <t  <  1 28  msec 

where  f1  =  230  Hz  f3  =  250  Hz  f5  =  270  Hz 
f2=  240  Hz  f4  =  260  Hz 
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Fig.  13-  Simulated  Signals  for  Passband  Analysis 
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Fig.  14A,  S/N  -  »  dB  { Frequency  Residue  ion  A!  ~  7r$  Hi: 


Fig.  14.  FFT  Absolute  Magnitude  of  Signal  (Fig.  13)  of  128  Input  Points  at  S/N  =  oo  dB, 
S/N  =8.47  dB,  and  S/N  =  3  96  dB  (Frequency  Resolution  of  Af  =  7.8  Hz) 
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Fig.  15.  Z-Transform  Absolute  Magnitude  Response  of  128  Input  Point  Simulated  Signals 
(Fig.  13)  Evaluated  Along  Contour  3  (Fig.  7)  for  101  Output  Points  (Passband  200 
through  300  Hz;  Frequency  Resolution  Af  =  2  Hz;  and  S/N  =  oo  dB 
S/N  =8.47  dB,  and  S/N  =  3.96  dB) 
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